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Abstract 



^ Recent mathematical fluid dynamics models have shed light into an outstanding 

^; problem in reproductive biology: why do spermatozoa cells show a 'preference' 

cn 

(_J for swimming near to surfaces? In this paper we review quantitative approaches 

'y—t 

^ to the problem, originating with the classic paper of Lord Rothschild in 1963. 

^ . , 

r^ A recent 'boundary integral/slender body theory' mathematical model for 

[ ^] the fluid dynamics is described, and we discuss how it gives insight into the 

, mechanisms that may be responsible for the surface accumulation behaviour. 

^ We use the simulation model to explore these mechanisms in more detail, and 

<^ . . . 

iy~\ discuss whether simplifled models can capture the behaviour of sperm cells. 



The far-fleld decay of the fluid flow around the cell is calculated, and compared 
with a stresslet model. Finally we present some new flndings showing how. 



O despite having a .ela.iveiy ™a„ i.yd.odyna™ drag, .he sperm cell 'head- has 

7— H very signiflcant effects on surface accumulation and trajectory. 
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1. Introduction 

Approximately one in six couples in the world experiences difficulties conceiving 
children, and on clinical examination male factors are found to contribute in at least 
half of all cases. A critical aspect of the function of the human male gamete is its 
motility: the sperm's ability to swim through the liquids of the female reproductive 
tract to the vicinity of the egg, penetrate the surrounding layers and fertilise. Motility 
is therefore an important target for novel diagnostic and therapeutic techniques, and 
advances in this area would have considerable financial and quality-of-life implications 
in human fertility, domestic animal breeding and wildlife conservation. The ability of 
sperm to swim is perhaps their most obvious function to the physical scientist, having 
long-provided an important archetype in mathematical fluid dynamics modelling. In 
this article, we discuss an important and surprisingly elusive aspect of the swimming 
behaviour of sperm both in vivo and in vitro — the phenomenon of 'surface accumula- 
tion'. 

Rothschild |24j famously observed the non-uniform, distribution of sperm in a drop 
of bull semen: spermatozoa typically exhibit a 'preference' for surfaces, with by far 
the greatest concentrations being found near to the inner faces of imaging chambers 
or the microscope slide and coverslip (Figure [l|. Moreover, once in the vicinity of a 
surface, sperm cells tend to remain at a constant distance, a fact that is very convenient 
in microscopical imaging of the sperm and its flagellar movement. Smith et al. |28j . 
In the sperm of humans and the majority of other species studied, gravity has not 
been observed to have a significant effect, so that the number of cells accumulating on 
upper and lower surfaces is typically similar, Winet et al. [29j. A related phenomenon 
occurs in other motile microscopic cells such as bacteria, and this has long interested 
researchers, particularly at the interface of mathematical modelling and experimental 
biology (see [2], [5], [E], [29]). 

A recent study of surface accumulation of sperm by Smith et al. [27] predicts that, 
at least for the human gamete, the effect can occur through fluid mechanics alone, 
resulting from the viscous-dominated regime of microscopic swimming, without the 
need for a specialised 'bias' of the flagellar beat to push the cell towards the surface. In 
this paper we review the findings of this study, then discuss how the model gives insight 
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Figure 1: Micrographs showing (a) high sperm cell concentration near a microscope slide, (b) 
low sperm cell concentration at a level mid-way between the microscope slide and coverslip. 
Unpublished data, Smith and Kirkman-Brown 2009. 

into the mechanisms that may and may not be responsible for this behaviour, and 
in particular some important and previously unexplored aspects that a fundamental 
theory of surface accumulation must address. 



2. Mathematical modelling of cell swimming at zero Reynolds number 

Fluid flow is modelled mathematically by the Navier-Stokes partial differential equa- 
tions. For a constant viscosity Newtonian incompressible fluid, these take the following 
form: 



p\ — + {u- V)m ) = - Vp + fiV^'u, V • w = 0. 



(1) 



The dependent variables are the velocity field u = u{x, t) representing the fluid velocity 
at location x and time t, and the pressure field p = p{x, t). The constants p and /i are 
density and viscosity, the latter quantity being measured in Pascal seconds (Pa s) in 
SI units. The first equation is derived from the principle of momentum balance, with 
the terms on the left hand side arising from inertia, while the terms on the right hand 
side arise from pressure and viscosity respectively. The second equation corresponds 
to the principle of mass conservation, combined with incompressibility. 

The most famous source of mathematical difficulty is the nonlinear term (u ■ V)tt, 
from which arises the great unsolved problem of turbulence and the open question of 
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whether the Navier-Stokes equations always have smooth solutions. However, in sperm 
fluid mechanics it is possible to neglect this term: nondimensionalising the equations 
with U as velocity scale, L as length scale, L/U as time scale and ^U / L/^ as pressure 
scale, the momentum equation is simplified to 

^(^ + («.V)t.)=-V^ + VV (2) 

where the coefficient on the left hand side is the Reynolds number R = pUL/ ^^ first 
identified by Osborne Reynolds as the parameter governing the transition to turbulence. 
A sperm cell has length approximately 50;um, swims with velocity approximately 
50 /im/s, and the density and viscosity of saline solution can be estimated as 10^ kg/m^ 
and 10"'^ Pa s. We estimate the Reynolds number as 2.5 x 10"'^, and indeed in solutions 
of higher viscosity it will be smaller still. It is therefore an excellent approximation to 
consider the 'zero Reynolds number' Stokes flow equations, 

= -Vp + V2w, V-M = 0. (3) 

Solving this equation directly using the computational fluid dynamics approaches 
generally applied to higher Reynolds number flows is computationally expensive, since 
there are three spatial dimensions. Simulating the trajectory of a sperm additionally 
involves solving the equation at a sequence of timesteps, and an additional complication 
is the presence of moving boundaries due to the movement of the flagellum and 
head, which in a computational fluid dynamics solution would necessitate moving the 
computational mesh at every timestep, which incurs further expense. The analysis we 
use in this paper exploits the linearity of equation ([3| to solve the problem much more 
efficiently. 

Very low Reynolds number corresponds to the near-absence of inertia, and in this 
respect, the fluid mechanics of sperm swimming differs significantly from more familiar 
situations such as human or fish swimming, or indeed most of the fluid dynamic 
phenomena we encounter. An important property of equation ([3]) is the absence of 
explicit time dependence, leading to time-reversal symmetry. Classic demonstrations 
of this were given by Sir Geoffrey Taylor in his 1967 fllm Low Reynolds number flow, 
including the 'smearing' and 'unsmearing' of an ink droplet in glycerine, and the 
abolition of swimming of a model fish by high viscosity liquid. This latter effect is 
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a crucial aspect of sperm swimming — only by performing a time-irreversible motion, 
such as propagating a bending wave, can sperm swim effectively. 'Flapping' motions of 
a fish tail or other reciprocating movements are ineffective — later referred to by Edward 
Purcell [23^ as the 'Scallop Theorem'. Recent research by Lauga 16J has shown that 
in non-Newtonian liquids, such as polymer suspensions, the introduction of an explicit 
time-dependence results in the breakdown of this theorem. 

3. Singular solutions of the Stokes flow equations 

The problem of sperm swimming will be formulated using equation ([3]), together 
with boundary conditions for the fluid velocity on the surface of the sperm tail, 
head, and additionally a nearby glass surface. The linearity of equation ([3| allows 
the superposition of solutions in order to satisfy the boundary conditions, and the 
techniques we shall use will be based on singular solutions, in particular the Stokeslet 
(point-force), Stokes-dipole (point-force dipole) and the potential flow source dipole. 

The Stokeslet Sjk{x,y) is defined as the solution u{x) = {Sik{x,y)., S2k{x,y), 
S3k{x, y)) and p{x) = Pk{x, y) of 

- Vp + V^w + efc(5(a; - y) = 0, V • w = 0, (4) 

which is the Stokes flow equation with a concentrated unit force per unit volume located 
at y, pointing in the ^-direction, with 5{x — y) denoting the Dirac delta distribution. 
The solution for a point force in an infinite fluid is given by Sjk = {Sjk/r + rjrk/r^)/8Tr 
and Pk = rj/{4:TTr^), where rj = Xj — yj and r^ = r\ + r2 + r|, and 5jk denotes the 
Kronecker delta tensor. Other singular solutions include the Stokes-dipole, defined as 
^fki — {d/dxi)Sjk, and the potential source-dipole PR = {d / dxk){rj /r^)) /Att . The 
latter is named due to being a source-dipole solution of the potential flow equations, 
and is a solution of equation ([3]) with zero pressure field. 

The Stokeslet, coined by Hancock |12j as a name for the purely viscous component of 
the Stokes flow around a translating sphere, the other component being the potential 
source-dipole. The Stokeslet and potential-dipole tensors are singular at the centre 
of the sphere, so that the flow is regular in the domain |a; — ^/l > o, where a is the 
sphere radius. The concept of representing flows around finite bodies by combinations 
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of singularities outside or on the boundary of the flow domain forms the basis for the 
slender body theory and boundary integral method we discuss below. 

4. The no-slip condition in viscous fluid mechanics and the method of 

images 

The appropriate velocity boundary condition in viscous flow is the no-slip, no- 
penetration condition u{x = X) = U{X), where X is the position vector of a solid 
surface and U{X) is its velocity. In this problem we are concerned with a stationary 
boundary representing the glass surface of a microscope slide or a coverslip, and a 
moving boundary representing the cell head and tail surfaces. The glass surface 
is modelled as the plane X3 = 0, with flow taking place in the 'half-space' region 
xs > 0. In the experiment, there is an additional bounding surface, but since we shall 
only consider motion close to one surface, this is neglected in our model. Liron and 
Mochon 20J showed how to model the effect of the additional surface, but this is more 
computationally expensive, as discussed in Smith et al. [5B1. 

Blake [3] addressed the problem of modelling a flow using Stokeslet distributions 
while satisfying the no-slip condition it (2:3 = 0), motivated by the propulsion of fluid 
by cilia protruding from a cell surface. The technique is equally useful in modelling 
the interaction of a glass surface and a sperm cell. We denote the solution here as B, 
where 
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(5) 



dRk 

with Ra = ra, R3 = 3^3 + ^3 ^^d the index a taking values 1, 2 only. Equation ([5]) 
consists of a Stokeslet in the fluid, an equal and opposite image Stokeslet, and further 
images consisting of a Stokes-dipole and potential dipole. The tensor {SjaSak — SjsSsk) 
takes the value -|-1 for j = k = 1,2, value —1 for j = k = 3 and zero otherwise. An 
equivalent solution had technically been found previously in the 1896 work of Lorentz 
(see [H]); however it was not until solution ([5| was discovered that the nature of the 
image systems was clear. Furthermore this result showed how the far-fleld is modifled 
by the boundary effect: the 0{r~^) decay of a force acting parallel to the boundary is 
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modified to become a 0{r~'^) decay, the far-field being of 'stresslet' character, i.e. a 
symmetric Stokes-dipole. This has recently been used to understand the fluid dynamics 
of the movement of vesicles in the developing embryonic node, Smith et al. [25j . 

5. Slender body theory and the boundary integral method 

Performing the first detailed mathematical analysis of microorganism swimming, G. 
J. Hancock |12j . working with Sir James Lighthill, formulated a slender body theory 
for Stokes flow, inspired by slender body theory for potential flow. This involved 
representing the fluid flow around the flagellum using a line integral of Stokeslets 
weighted by unknown force density. Considerable work on simpliflcations to slender 
body theory (see [TU], HE]); adaptations for cilia modelling (see [i], 0, [H], P^). 
refinements to incorporate a spherical cell body (see ^^), large-amplitude motions 
(see [H], [13]), asymptotic analysis (see [I], [0]) and higher-order corrections (see [T3]) 
have been presented, as reviewed in detail in by Smith et al. [57] • Essentially, all of 
these theories involve representing the fiow field due to the moving tail by a line integral 



Jo 



tttaii = / G{x,^) ■ f{s) ds, (6) 

where the kernel G is a force singularity Green's function, coupled with a potential- 
dipole P to improve the accuracy of the slender body representation, and f{s) is the 
force per unit length. 

In order to model the effect of the sperm cell body or 'head', there are three fluid 
dynamic effects that must be taken into account. In order of importance, they are: 

1. The drag and torque on the translating and rotating head, balancing the force 
and torque produced by the flagellum. 



2. The near and far-fields of the fiow generated by the head as it moves through the 
fluid, 

3. The fact that the no-slip condition must be satisfled on the head, taking into 
account the flow produced by the flagellum. This modifies f{s) and usually 
reduces the swimming speed. 
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The only fluid dynamic model to have taken all these effects into account while using 
slender body theory exclusively was that of Higdon |13j , which modelled the swimming 
of a spherical-headed microorganism. Taking into account the non-spherical head of 
most sperm cells, and additionally the effect of a nearby no-slip surface, requires the use 
of more complex numerical implementations, such as the boundary integral method, as 
first applied to sperm by Phan-Thien et al. 22J. The boundary integral method uses 
distributions of Stokeslets over the a surface dH, with force per unit area given by the 
unknown traction (f){X , t) for X G dH. Recently, we used a hybrid of the boundary 
integral method for accurate modelling of the head, and slender body theory for the 
flagellum (see [17]), in order to provide an efficient and accurate representation of both 
parts of the cell. In this method, the velocity field for points in the fluid x outside of, 
or on, the surface of the sperm is given by 

u{x,t) = Mhead(aJ,t) + iitaii(a;,t) 

= // B{x,X'^).cl){X\t)dX'^+ f G{x,^is,t)).f{s,t)ds. (7) 

J JdH Jo 

Since the movement of the flagellum will be specifled relative to the cell, not a sta- 
tionary observer, a coordinate transformation is required. The position of the flagellar 
centreline in the laboratory frame is given by 

^ = Xo + M-^', (8) 

where Xq is the location of the head/flagellum junction, M is a rotation matrix from 
the 'body frame', in which the head is stationary, to a flxed 'lab frame', and ^ is the 
flagellum position vector in the body frame. The flagellar movement is specifled as the 
function ^ . The body frame will translate and rotate with a priori unknown velocity 
V and angular velocity 17. By rigid body mechanics, with respect to the laboratory, the 
kinematic velocities of points on the head U{X ) and flagellum U{^) are respectively 

U{X'') = V + rtA{X^ -Xo), 

C/(|) = V + nA(|-Xo) + M.^', (9) 

In order to close the system so that (j), /, V and Jl can be determined simultaneously, 



Surface accumulation of spermatozoa 9 

we use the principles of force and torque balance in zero Reynolds number flow: 

= [ [ 4>iX\t)dX''+ f fis,t)ds, 

J JdH Jo 

(X'^ - Xo) A ct,{X\ t) dX'^ + / (|(s, t) - Xo) A f{s, t) ds. 
an Jo 

(10) 



At each instant in time, the initial position and orientation are known, and the 
fiagellum position ^, and kinematics ^, with respect to the head are given by a mathe- 
matical formula chosen to represent the flagellum waveform. From this, equations ([7]) 



and ( 10 ) form a closed system for the unknowns (p, /, V and fi. A collocation method 
is used to discretise these equations spatially, and a matrix system with typically 
3 X (60 + 32) + 3 + 3 scalar unknowns is formulated, corresponding to 60 constant-force 
segments for the fiagellum, 32 constant-force elements for the head, the additional 
terms being the 3 scalar components of V and 17. The solution of this system allows 
the updating of the cell position and orientation from V and 17. From the new values 
of / and cf), the equation ([t]) can be used to calculate the velocity of the fluid at any 
point in space. 

This method is very computationally efficient, allowing the accurate simulation of 
8000 beat cycles of migration, discretised with 100 timesteps per beat, in approxi- 
mately 20 hours of CPU time. Such efficient computation of a moving-boundary, 3D, 
unbounded and time-dependent problem is not possible with standard computational 
fluid dynamics methods that involve the determination of the fluid velocity at every 
point in the domain at every timestep. The simpliflcation arises from the linearity of 
the Stokes flow equations, the corresponding ability to take integral sums of known 
solutions, and our knowledge of Green's functions such as B. 

Simulations were performed with the University of Birmingham's BlueBear cluster, 
using a single core per simulation. This allows parameter space to be searched through 
the simultaneous execution of many simulations; however a single simulation could be 
performed in a similar time period on a desktop PC. 
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6. Prescribing the flagellar beat 

The remaining step in the model is to specify the movement of the fiagellum £ . 
While some approaches have generated the flagellar beat from a complex interaction 
of internal mechanical activity and external fluid forces, we have taken the simpler 
approach of writing ^ as an analytic function. The flagellar waveform of human sperm 
is usually a three-dimensional shape, likely to resemble a nearly-planar flattened helix 
(see [14] ) . In [27] we examined the effect of both a planar waveform and a flattened 
helix, and found that the difference in surface accumulation in the two cases was 
minimal. Hence in what follows we shall focus on the planar beat, which simplifies the 
analysis of the fiow field around the cell. The waveform we choose was first suggested by 
Dresdner and Katz i7|, and is proportional to a sine- wave sm{kxi—ujt), modified to have 
linearly increasing amplitude along the fiagellum. In [27J , we found that the angular 
wavenumber parameter k is the critical parameter governing surface accumulation, as 
discussed in the next section. Figure [2] shows the waveform for k = 2tt and A; = Svr 
respectively. 

(a), k = 2tt (b), k = Svr 




Figure 2: Planar sinusoidal waveforms used for the simulations, with (a) fc = 27r and (b) 
k = Stt. 



7. Recent findings 

Our recent simulation study [27] mainly concerned the behaviour of a cell initially 
at X3 = 1 in non-dimensional units, i.e. one body length away from a no-slip boundary 
at X3 = 0, with initially parallel orientation, and its beat plane parallel to the no-slip 
boundary, as shown in the simulation snapshot Figure Isj^a). The behaviour predicted 
by the simulations for wavenumbers k > 2 Air is shown schematically in Figure [31 
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the cell initially pitches towards the surface and swims towards it, then pitches away, 
repeating this behaviour until it approaches an 'equilibrium height' at which it remains. 
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Figure 3: (a) The initial configuration for all simulations, with the head/flagellum junction 
initially at a; = (0, 0, 1). (b) Schematic showing the surface accumulation behaviour predicted 
by simulations for wavenumbers k > 2.47r. The horizontal distance travelled is an order of 
magnitude greater than shown in the schematic, as shown in Figure^ 

Figure llFa) shows the definition of height, measured from the surface to the head/ 
flagellum junction, and pitch angle, measured so that a positive value corresponds to 
pitching away from the surface. Simulation results for height and pitch angle with 
wavenumbers k = 2tt and fc = Svr are shown in Figure [4J These results are based 
on those presented in [27], which were the first long timescale simulations of sperm 
behaviour near surfaces, although the results we present here are for a greater number 
of timesteps. 

These results show that accumulation can occur purely through fluid dynamic 
effects: Figure He) shows a relatively rapid convergence of a cell with wavenumber 
A; = Svr to a specific height above the surface — around X3 = 0.31. Figure [4|b) shows 



12 D. J. SMITH AND J. R. BLAKE 

the behaviour of a ceh with wavenumber k = 2tt. While it does not converge, it does 
suggest that after an initial period of swimming away, the cell will eventually return, 
swimming towards the surface. It is unclear whether this could ever result in swimming 
at a stable height, as for k = Svr, and indeed whether the oscillations will grow or decay 
in amplitude. A simulation period of 8000 timesteps corresponds to an experimental 
timescale of 400-800 seconds, depending on the beat frequency. 
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Figure 4: (a) Schematic showing the definitions of height above the surface, measured to the 
head/flagellum junction, and pitch angle, measured from a plane parallel to the surface, (b) 
Height and pitch angle trajectories, from a simulation with wavenumber k = 27r. (c) Height 
and pitch angle trajectories from a simulation with wavenumber fc = Stt. 



No specialised 'bias' in the flagellar beat towards the surface is required to produce 
the accumulation behaviour, provided the wavenumber is sufficiently large. As observed 
experimentally by Winet et al. |29J , cells eventually swim at a finite constant distance 
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from the surface, as opposed to necessarily swimming against it. The surface does 
not simply attract the cell; it causes alternate pitching towards and away from the 
surface that steers the cell to this finite distance. Finally, perhaps counterintuitively, 
the 'equilibrium state' of a cell swimming parallel to a surface involves a slight angle 
of inclination away from the surface. In the following sections we investigate further 
the fluid mechanics of this behaviour. 



(a) 



(b) 




no-slip surface 



(c) 
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j 'attraction field' 



resultant trajectory 



Figure 5: (a) Schematic showing the planes parallel to the no-slip surface used for the 
calculation of the mean flow field around the cell.(b) Density plot showing the vertical 
component of the flow field us, averaged over one beat cycle, and averaged between the 
planes Zo — s and Zo + e, where e = 0.02. (c) The equilibrium parallel-swimming behaviour 
of the cell involves a balance of inclined swimming away from the surface, and the attraction 
flow field shown in (b). 

8. The attraction field and pitching behaviour 

In order to gain insight into this effect, we calculate the time- averaged flow field 
around the cell, in particular the vertical component u^ that may push the cell towards 
or away from the surface. The flow fields evaluated at any finite distance above and 
below the cell show almost equal and opposite behaviour, so instead we calculate the 
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average ^ {u{Xo + £63) + u{Xo — £63)), as shown in Figure pTa). In an infinite fluid, 
this average would be zero, however the image singularities of B (equation (pi)) generate 
a non-zero vertical component, the dominant image singularity being the Stokes-dipole. 

An example result for Xq = (0,0, 1) and A; = Svr is shown in Figure pTb), showing 
that the cell generates an 'attraction held', i.e. a vertical fluid velocity that pushes it 
towards the surface. This occurs for wavenumbers k = 2tt and A; = Svr, and also occurs 
very close to (X3 = 0.1) and further from [xj, = 1.0) the surface. In accumulating 
cells, the eventual equilibrium trajectory occurs due to a balance of the attraction 
field, drawing the cell towards the surface, and a progressive component away from 
the surface, occurring because there is a small pitch angle away from the surface — see 
for example Figure ETc), where the eventual pitch angle converges to 0.42° while the 
height converges to X3 = 0.31. This mechanism is summarised in Figure iMc). 

9. The propulsive force distribution on the flagellum 

Cisneros et al. [5], using the related method of regularised Stokeslets, modelled the 
swimming of a bacterium near to a wall, representing the head as a sphere and the 
flagellum as a propulsive stick, with constant propulsive force distributed along its 
length. Surface accumulation behaviour was interpreted from the flow fleld generated 
around the swimmer, and also the existence of a resultant force predicted to act on 
a swimmer constrained to move parallel to a surface. In our simulations of sperm 
swimming we take a different approach: the propulsive force is calculated directly 
from the fluid/flagellum interaction, and cells are predicted to swim parallel to the 
surface while remaining free from any resultant force. 

While the attraction field appears to occur regardless of the wavenumber and dis- 
tance from the surface, the pitching behaviour is very sensitive to both of these 
parameters, and it is this pitching effect, rather than the attraction field, that is 
responsible for the behaviour shown in Figure p[b-c). Wavenumber affects the flow 
field through changes to the force distribution on the fiagellum. In order to quantify 
how this changes with wavenumber, we calculate the time-averaged force distribution 
resolved in the direction of propulsion. Graphs of this quantity for k = 271 and A; = Svr 
are shown in Figure [6J both of the time-averaged force distributions differ significantly 
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from the constant function used in the bacterial model. The average propulsive force 
along the flagellum, which balances the drag on the head and causes progression, is a 
fraction of the peak magnitude of the force along the flagellum: for a swimming cell 
such as a sperm with a relatively small cell body, the vast majority of drag occurs 
along certain portions of the flagellum, rather than the head. 



(a), k = 27r 



(b), k = 37r 
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Figure 6: Flagellar force distribution, averaged over one beat cycle for (a) fc = 27r and (b) 
k = Stt. The dotted line denotes the mean force distribution along the flagellum, showing 
that most of the propulsive effect generated by the mid region of the flagellum is balanced by 
nearly-equal and opposite drag occurring at the ends, as opposed to the head of the cell. 



A fundamental theory of surface accumulation must take into account the complexity 
of the flagellar force distribution. Indeed, the relatively small change between k = 2tt 
and fc = Svr results in very significantly different behaviour (Figure l4| . It should be 
noted that the bacterial flagellum differs substantially in structure and waveform from 
the sperm flagellum, and the force distribution for this cell type is likely to differ from 
our results in Figure [6] — and may more closely resemble a constant function. 

The simplest model of the fluid flow produced by a force- free swimmer is a stresslet, 
i.e. a Stokes-dipole that approximates a propulsive force due to the flagellum and a 
drag force due to the head acting very close together. This is based on the following 
simplified model: the drag on the head pushes fluid in front of the cell in the direction 
of swimming, while the propulsive effect of the flagellum pushes fluid rear of the 
cell in the opposite direction. Analysis of the propulsive force distribution Figure [6] 
shows that the situation is more complex for sperm motility. Figure [6] leads us to 
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examine the hypothesis that the flow may be modefled more accurately as a Stokes- 
quadrupole {d'^ /dx\)Sji, at least for fc = Svr, as shown in Figures] This singularity 
has far-field decay 0{r~^), compared with the 0{r~'^) behaviour of the Stokes-dipole. 
From our simulations, the decay of the velocity field around a sperm in an infinite 
fluid is approximately 0{r~'^-'^), intermediate between the Stokes-dipole and Stokes- 
quadrupole. When the no-slip surface is present, the decay of the velocity field is even 
more rapid, being approximately 0(r~^'*), as may be expected from presence of the 
additional image systems. 

(a) (b) 




drag propulsion drag ^^ ^^ 
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Figure 7: Two simplified models for the force exerted on the fluid by a sperm: (a) Stokes- 
dipole, (b) Stokes-quadrupole, the latter being inspired by the force distribution shown in 
Figure l6[c). 

10. The effect of the cell 'head' 

Figure [6] suggests that the drag on the head, which is equal and opposite to the 
mean of the propulsive force on the flagellum, is very small in comparison with the 
peak flagellar force. This leads us to test whether the head has any significant effect on 
the accumulation behaviour. We perform simulations of a cell with half the usual head 
dimensions, resulting in l/8th the volume, and also a 'pinhead' cell — as sometimes 
observed experimentally — with volume 1/lOOOth that of the normal head. The results 
are shown in Figure [S] — a reduction in head dimensions results in much slower decay of 
the oscillations, while a pinhead cell does not even approach an equilibrium trajectory 
over 8000 beat cycles. Despite the fact that the head makes a relatively minor contri- 
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bution to the drag in comparison with regions of the flagenum, it is nevertheless crucial 
to the surface accumulation effect. The reasons for this are unclear, but emphasise the 
subtlety of the fluid dynamic interaction. 
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Figure 8: Simulation results showing the effect of the head proportions on accumulation 
behaviour, (a) A cell with head dimensions half that of a normal cell, giving l/8th the 
volume. This cell is predicted to accumulate, but the oscillations decay much more slowly 
than in the normal case, (b) A 'pinhead' cell with head dimensions one tenth that of a normal 
cell, giving 1/lOOOth the volume. This cell is not predicted to accumulate over the timescale 
simulated of 4000 beat cycles, but nevertheless does exhibit an eventual pitching towards the 
surface. 
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11. Conclusions 

Despite the apparent simplicity of the problem — modelling a single cell swimming in 
zero Reynolds number flow — surface accumulation is a relatively subtle problem and 
has evaded a simple explanation since the 1960s. Fluid dynamic simulations give a 
number of insights into the mechanism for this behaviour, and the roles of the head 
and flagellum force distributions. To summarise: 



For sufficiently large wavenumber, sperm have a 'preferred' trajectory, in which 
they swim near and parallel to a surface, with zero total force. 

This occurs through a combination of pitching behaviour, which 'steers' the cell 
to an equilibrium height and a relatively weak attraction field, that balances the 
eventual slight inclination of the cell away from the surface. 

Any theory of sperm accumulation must take into account the shape of the 
propulsive force distribution, and the fact that the sperm does not behave as 
a simple Stokes-dipole drag/propulsion combination. 

The sperm head, while contributing a relatively small fraction of the drag, has 
very signficant effects on accumulation, in particular causing more rapid conver- 
gence of the cell trajectory. 

In the real physiological system, sperm undergo transient spontaneous changes 
in swimming behaviour. This may be responsible for the fact that observations 
(for example. Figure 1) rarely show every single cell swimming stably at a 
fixed distance from the surface, as might be predicted from simulation studies. 
Furthermore this is an important difference between sperm and other swimming 
cells such as bacteria. 



A 'simple' explanation of this long-standing problem is surprisingly elusive. Moreover, 
similarities and differences between sperm and bacterial behaviour near surfaces require 
further investigation. 
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